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The cross-correlation of a foreground density field with two different background convergence fields 
can be used to measure cosmographic distance ratios and constrain dark energy parameters. We 
investigate the possibility of performing such measurements using a combination of optical galaxy 
surveys and HI intensity mapping surveys, with emphasis on the performance of the planned Square 
Kilometre Array (SKA). Using HI intensity mapping to probe the foreground density tracer field 
and/or the background source fields has the advantage of excellent redshift resolution and a longer 
lever arm achieved by using the lensing signal from high redshift background sources. Our results 
show that, for our best SKA-optical configuration of surveys, a constant equation of state for dark 
energy can be constrained to ~ 8% for a sky coverage faky = 0.5 and assuming a (T(f2DE) = 0.03 
prior for the dark energy density parameter. We also show that using the CMB as the second source 
plane is not competitive, even when considering a COrE-like satellite. 


I. INTRODUCTION 

Cross-correlation cosmography is a method that uses 
the gravitational lensing of galaxies in different source 
(background) redshifts in correlation with a given lens 
(foreground) population in order to constrain dark energy 
W- If the foreground distribution is sufficiently narrow 
in redshift, the ratio of the foreground galaxy density (5g) 
- background lensing convergence (k) cross-correlation 
from different background bins is a purely geometrical 
quantity which depends only on the source distribution 
kernels. In the case that the background distributions 
are also narrow in redshift, this quantity simplifies to a 
ratio of comoving radial distances. 

The distance ratio depends on the background cosmol¬ 
ogy, and in particular on the amount of dark energy and 
its equation of state; it is independent of the large scale 
structure details (power spectrum, bias) and the angu¬ 
lar scale. Therefore, measurements at different angular 
scales - even those in the non-linear regime - can be used 
to constrain dark energy. Of course, removing the de¬ 
pendence on the growth rate and the associated informa¬ 
tion it carries about dark energy means that the derived 
constraints will be weaker than the ones derived via the 
usual methods which depend on both the geometry and 
the growth, such as lensing tomography [SHE]. However 
the geometric method directly probes the evolution of 
the Universe via the redshift-distance relation and it is 
therefore a robust consistency check. Furthermore, iso¬ 
lating geometry from structure formation is very useful 
for General Relativity tests, since modified gravity mod¬ 
els can mimic the expansion history of the Universe but 
have different clustering properties [^. 

In this work we investigate these correlations using a 
combination of HI intensity mapping surveys and galaxy 
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surveys. Intensity mapping is an innovative technique 
able to map the large-scale structure of the Universe in 
3D (see e.g. [10] )• K uses HI as a tracer of the dark 
matter density field by measuring the intensity of the 
redshifted 21cm line over the sky in a range of redshifts 
without detecting individual galaxies, treating the 21cm 
sky as a diffuse background, similar to the CMB. In 
m the potential of the planned Square Kilometre Array 
(SKA) [T^ to deliver HI intensity mapping maps over a 
broad range of frequencies and a substantial fraction of 
the sky was investigated. Using HI intensity mapping 
to probe the foreground density tracer field and/or the 
background source fields has the advantage of excellent 
redshift resolution and a longer lever arm achieved by us¬ 
ing precise measurements of the lensing signal from high 
redshift background sources. 


The paper is organised as follows: In Section |TT] we de¬ 
scribe the general formalism we will use. We form the 
tracer density-convergence power spectrum correlating 
the density from a foreground bin (/) and the conver¬ 
gence from a background bin (b). We then derive the 
distance ratio, describe how it depends on dark energy 
and formulate the corresponding Fisher matrix for the 
derivation of dark energy constraints. We then move on 
to the results. Sect ion [iHj We first consider the usual case 
of a galaxy survey for both the foreground lenses and the 
background sources, employing Large Synoptic Survey 
Telescope (LSST)-like parameters [T3]. Then we consider 
combinations of optical galaxy surveys (LSST) with HI 
intensity mapping surveys performed with SKA-like in¬ 
struments for measuring the foreground density fields and 
the background convergence fields. Finally, we show how 
an even longer lever arm could be achieved by considering 
21cm lensing from the Epoch of Reionization, or CMB 
lensing with a future satellite mission. For all cases con¬ 
sidered we derive constraints on the constant equation of 
state parameter wq and the dark energy abundance Hde- 
We discuss our results and conclude in Section HVl 
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II. FORMALISM 


We are interested in the correlation of the foreground 
(/) density tracer field <5tr with the background (5) con¬ 
vergence K field. We will assume linear biasing, so that 
i^tr = bti-S, where 6 is the underlying matter density held 
and btT is the tracer’s bias — in this work, we use ‘tr = 
for galaxy surveys and ‘tr = HI’ for intensity mapping 
surveys. Using the Limber approximation |14j the angu¬ 
lar cross-power spectrum Cs^^k, is given by 

CsUe; f,b) = J ^^Wfixf) JdxbW.ixb) 


where x is the comoving distance and Pss is the three- 
dimensional dark matter power spectrum. We have as¬ 
sumed that the foreground galaxy distribution Wf and 
the background galaxy distribution Wb do not overlap, 
so that Xf < Xb always.This assumption is important, as 
a possible contamination of background sources in front 
of foreground density tracers can be a problem; the fore¬ 
ground populations have to be completely in front of the 
background populations for Eq. Q to be exact. A com¬ 
prehensive analysis of this effect for the case where both 
the foreground and background populations are probed 
by galaxy surveys is presented in Of course, one 
can avoid this systematic by selecting the foreground and 
background populations such that they are clearly sepa¬ 
rated in redshift. Using HI intensity mapping is advanta¬ 
geous, because it provides very precise redshift informa¬ 
tion which is not the case when using photometric galaxy 
surveys. 

If the foreground lens slice is narrow enough in red¬ 
shift {Az ~ 0.1 is sufficient), we can approximate the 
foreground redshift distribution as a delta function at a 
distance x/, Wf{xf) = ~ Xf)- We then find 


f,b) 


l btAXf)Pss[^,Xf) 
2c2 I a(Xf)Xf 


X f dxbWbixb) X (2) 

J Xb 


By taking the ratio of measured using two different 
background distributions bi , 62 but the same foreground 
/, we see that the details of the large scale structure and 
the bias cancel out and all that remains is the geometrical 
ratio of the source distribution kernels. In the limit that 
the two background sources also have a delta-function 
distribution (at Xbi, Xb 2 > respectively) the ratio becomes 


R{zf,Zb,,Zb^) 


/,&i) 

f,b2) 


(Xbl - Xf)/Xbl 

(Xb 2 -X/)/Xb 2 ' 


Note that the ^-dependence (angular scale dependence) 
also cancels out. Assuming a flat Universe, the right- 
hand-side depends only on the dark energy parameters 


{fluE,wo,Wa), where we have used the parametrisation 
w = Wo + Wa{i — a) [15] for the dark energy equation of 
state. 

The ratio R defined in Eq. ([^ is observable and de¬ 
pends only on the comoving distances to the foreground 
and background populations. In a flat Universe, the co¬ 
moving distance is given by 

x(z) = ^ rdz'[(l-HDE)(l + ^')' 

Jo 

+ Hde (1 + ^') 3 (i+“’t.+’"Ue- 3 >^a^ 7 (i+^')]-i/ 2 ^ ( 4 ) 


We shall now derive dark energy constraints using the 
Eisher matrix 


F, 


pq 


I 

pEE A^/sky(2£+l) 

/ ■^min 


dR^Idp 

Rf 


1 

W) 


dRf/dq 
Rf ’ 


where Rf = R{zf, Zb^j Zb^) is the ratio corresponding to 
foreground bin /, {p,q) = (Hde,' ib'Oj w^a) and o'^(I') is the 
fractional variance of Rf given by [3| 


+ J^bM] + [CsUfiJ,b2)]-^[Cs,^sUf,f) +A//W)] 
[a.(7 b2) + MbM)] - /, /, b2)-^ 

FMf{i)\CU^M,b2)- ( 6 ) 


Here, J^f{t) is the noise in the foreground density power 
spectrum f) measurement, and Nbi{(.) is the 

noise in the convergence power spectrum C'kk(I’, bi) mea¬ 
surement. The exact form of these noise terms depend 
on the kind of survey we are using (see next Section). 
Einally, /sky is the fraction of the sky our chosen surveys 
map, and A£ the binning in Aspace. 

The above equation is the Fourier space equivalent 
of the variance used in [ 1 | where halo-shear correlations 
were considered, but in that work only the shot variance 
contribution was used, i.e. the approximation 


a}{e) :^[CsU^J,bi)]-^Mf{e)MbAe) + [Cs,^.{i,f,b2)]-^ 

(7) 

was made |3] . Throughout this work we use the full noise 
expression of Eq. Q; however, we will ignore the (weak) 
covariance between different foreground bins (see |3] for 
details). 

Before we derive our results using different combina¬ 
tions of intensity mapping and galaxy surveys, it is useful 
to demonstrate the dependence of the ratio R on the dark 
energy equation of state. In Figurej^we show the depen¬ 
dence of R given by Eq. ^ with Zf = 0.5, Zb^ = 1.3 
and Zb 2 = 2.4 on the dark energy parameter wo (assum¬ 
ing Wa = 0) with Hde = 0.73. We immediately see that 
the dependence is very weak, R oc which means 

that in order to have good constraints on wo we would 
need a very accurate measurement of R. 
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FIG. 1: The ratio R as a function of wo for Zf = 0.5, = 1.3 

and Zb 2 = 2.4, assuming Wa = 0 and Ude = 0.73. The depen¬ 
dence on the equation of state is very weak {R oc 


III. RESULTS 


In the following we will investigate a number of dif¬ 
ferent survey combinations to perform cross-correlation 
cosmography studies. Initially we consider a galaxy sur¬ 
vey with LSST-like parameters for both the foreground 
lenses and the background sources populations. Then we 
consider intensity mapping surveys using SKA-like in¬ 
struments for measuring the foreground density field or 
the background convergence field, and compare them to 
the LSST-only results. We also consider the case where 
our second ( 62 ) background sources are at the Epoch of 
Reionization. Finally, we use the CMB as the second 
source plane and a COrE-like satellite [16]. 

When using Eq. ([^ for our Fisher matrix calculations 
we apply the delta function approximation for the dis¬ 
tribution of the foreground populations, as they are thin 
enough in redshift (Az < 0.1) in all cases. In the case 
where the background sources are probed with an opti¬ 
cal galaxy survey, we use the known redshift distribution 
of these galaxies. In the 21cm lensing case using the in¬ 
tensity mapping method, the HI galaxies distribution is 
unknown. However, our chosen bands are narrow enough 
so that the delta function approximation is valid (i.e. an¬ 
gular distances vary very little within each band). The 
delta function approximation is also valid when CMB is 
used as the second source plane. 

The survey combinations we considered together with 
the best obtained Icr uncertainties on wq for each case 
are shown in Table HI 


Section/Surveys 

/ 

61 

^2 

Swo 

ihIa/b 

LSST 

LSST 

LSST 

0.12/0.10 


III 

B 

SKA-like 

LSST 

LSST 

0.08 


III 

c 

LSST 

LSST 

SKA-like 

0.08 


III 

D 

LSST 

LSST 

SKA-like|EoR 

0.11 


III 

E 

LSST 

LSST 

COrE-like|cmb 

0.22 


TABLE I: The survey combinations considered in this work. 
For each combination the best obtained Icr marginalised error 
on the constant dark energy equation of state wo is shown. 


A. Results using a LSST-like survey for the 
foreground lenses (/) and background sources (&i,62) 

In the case of galaxy surveys we have Str —>■ Sg and, 
consequently, Cs^k for the cross-correlation power spec¬ 
trum, with Sg = bgS and bg the galaxy bias. The noise 
terms are the shot noise contributions given by 

Af/(«) = 4; = (8) 

Uf 

where h® is the number density of foreground galaxies in 
bin /, is the shape noise of each background galaxy 
and hf {i = 1 , 2 ) is the number density of background 
galaxies in source bin i. 

We consider an LSST-like galaxy survey m with 
redshift range 0 < z < 3, galaxy number density 
40 arcmin“^, cr„ = 0.3, covering half of the sky /sky = 0.5. 
The redshift distribution of galaxies is modelled by 

^ oc (9) 

dz 

with zo = 1. We will consider 10 foreground top-hat 
redshift slices from 0.1 < z/ < 1 (i.e. Az = 0.1) and two 
source (background) populations between z = 1.1 and 
1.5, and between z = 2.2 and 2.6 - note that the first 
background source bin ( 61 ) will always be chosen to be 
close to the foreground population, as this leads to higher 
sensitivity of the ratio to the dark energy parameters, 
hence better constraints. Various auto, cross and noise 
power spectra that enter the calculation of the error in 
Rf given by Eq. ^ are shown in Figure using the 
Zf = 0.5 foreground (lens) bin and our two background 
(source) bins with Zf,j = 1.3 and z^^ = 2.4. In the Fisher 
matrix we take bins of A£ = 30 starting from = 45 
and going up to fmax = 3000. To model the various 
power spectra that enter the error calculation we use the 
non-linear fitting formula by Smith et al. HZI and we set 
the galaxy bias bg = 1. 
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FIG. 2: Various power spectra that contribute to the error in 
the measurement of , using the Zf = 0.5 foreground (lens) 
bin and two background (source) bins with = 1.3 and 
Zb2 = 2.4. 


Following previous works (see, for example, mm) , we 
will present constraints on a constant equation of state 
w = Wq employing a prior (T(nDE) = 0.03 for the dark 
energy density parameter. For the LSST-like survey we 
find 


6wo ~ 0.12 



^DE 

FIG. 3: Dark energy Icr and 2a constraints on wo and Dde 
fixing Wa = 0 and employing a prior on Dde correspond¬ 
ing to (j(Ddb) = 0.03. The fiducial model has w = —1 and 
Dde = 0.73. We have used an LSST-like galaxy survey with 
parameters described in the main text. 


As expected, the constraints we found are weaker than 
those from methods using the information encoded in the 
growth function. But the advantage of the geometrical 
method is that it carries less theoretical assumptions and 
can be used as a consistency check against methods that 
use both geometry and growth to constrain dark energy. 
For a detailed discussion on the growth-geometry con¬ 
straints and issues like photometric redshift accuracy and 
the use of a more general source scaling of lensing signals 
which works for both galaxy-shear and shear-shear cor¬ 
relations, see [ 3 ]. However, note that the aforementioned 
scaling relation is disrupted by the fact that we can only 
measure reduced shear and not shear, which introduces 
a scale dependence on the ratio making it somewhat sen¬ 
sitive to the matter power spectrum m- 

We will now move on to the interesting possibility 
of using intensity mapping surveys in combination with 
galaxy surveys to perform cosmographic studies. We con¬ 
centrate on the combination of LSST and SKA-like in¬ 
tensity mapping surveys. Note that possible synergies 
between the LSST and SKA for cosmological probes like 
galaxy clustering, weak lensing and strong lensing have 
been recently discussed in m- 


for the marginalised Icr wq uncertainty, while the joint 
(Hdej^o) Icr and 2cr contours are shown in Figure]^ 
Note that if we use the approximation for the squared 
fractional variance of R used in n, i.e. Eq. Q , we find a 
much smaller Sw ~ 0.04. This is in qualitative agreement 
with the findings in [3] and shows that, even if there is 
partial cancellation among the extra terms in the full 
formula presented in Eq. (§, they do not cancel exactly 
and should be taken into account. 


B. Results using a SKA-Mid-like Intensity 
Mapping survey for the foreground lenses (/) and a 
LSST-like survey for the background sources (61,62) 

In the case of intensity mapping (IM) surveys we have 
<)tr —>■ <)iM with (5 im = TSm = Tbm^, where T is the mean 
brightness temperature and 6 hi the HI bias. In a recent 
publication [20] the majority of current and planned HI 
intensity mapping experiments at redshifts z < 4 were 
analysed. Intensity mapping at radio frequencies has a 
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number of attractive features, a very important one being 
the fact that we can automatically measure redshifts with 
extremely high precision, which circumvents the photo-z 
uncertainties that can be a serious source of systematic 
error in cross-correlation cosmography. 

Focusing on the SKA-Mid instrument, we can consider 
an intensity mapping survey in which the instrument op¬ 
erates in either “single-dish autocorrelation” or “inter¬ 
ferometer” mode [20]. We will use the latter for our 
analysis since the single-dish mode does not probe the 
small scales due to the limited angular resolution which 
is fixed by the size of the dish. On the other hand, the 
angular resolution of the interferometer is fixed by the 
array’s maximum baseline, hence it can probe very small 
angular scales (equivalently, high multipoles t). 

In this case, the noise contribution is dominated by 
the thermal noise of the instrument which is calculated 
using the formula (note we use the uniform distribution 
approximation) [2TH23| 




N 


A/tobs/cover^max(j^)^ ’ 


( 10 ) 


where Tgys is the system temperature. A/ is the chosen 
frequency (equivalently, redshift) window, and fobs is the 
total observation time; f'max(j^) is the highest multipole 
that can be measured by the array at frequency v (wave¬ 
length A), and is related to Dteh the maximum baseline 
of the core array, by ^max(A) = 27rAei/A. /cover is the 
total collecting area of the core array, Acou divided by 

^(Ael/ 2 ) 2 . 

For the foreground lenses (/) we will consider a SKAl- 
Mid-like configuration (i.e. SKA-Mid Phase 1). For the 
background sources ( 6 i,& 2 ) we consider the LSST-like 
survey that was described in the previous Section. We 
split the foreground population in bins of A/ = 40 MHz, 
covering the redshift range 0.1 < z < 1. The rest 
of the parameters are Acou = 0.08km^, Aei = 8 km, 
= 20 + 66 X ( 300 ^)”" "" K, and Us = lO^hrs 
[24| . The observed IM auto power spectrum is given by 


(^. /) = [T(z/)]2&|iU5 (£, /) + QN, (11) 


where we used the fact that the redshift slices are thin 
and the mean brightness temperature 


r(z) = 566 


{l + zY 


0.003 J H{z)/Ho 


HK 


( 12 ) 


varies slowly within each bin. Therefore Nf(£) is given 
by (setting 6 hi = 1 ) 


Nf{e) = 


/^N 


[TizfW 


(13) 


We first compare the performance of combining SKAl- 
Mid (for the foreground lenses) and LSST (for the back¬ 
ground sources) with the performance of the LSST-only 
case presented in the last subsection. In order to model 


the various power spectra we use the non-linear fitting 
formula by Smith et al. m and set the HI bias 6 hi = 1; 
we also assume a constant Hhi = 4.9 x 10““^ which is 
the value calculated from the local HI mass function 
measured by HIPASS |2S|. For the same assumptions 
of the LSST analysis (in particular Uax = 3000), re¬ 
placing the LSST low redshift density measurement with 
SKAl-Mid observations provides constraints very consis¬ 
tent with those seen in Figure 

In order to utilise the smaller angular scales probed we 
next push the analysis to a larger value of f^ax to 10“^. As 
we can see in Fig. [^ the results are improved. We also 
present results considering the SKA2-Mid case (i.e. 10 
times bigger collecting area than SKAl) and £max = 10"^. 
At very high i values the shape noise terms dominate, 
and taking £niax > 10 "^ does not significantly improve the 
constraints. We should also note that since we use power 
spectra modelling to calculate the variance terms, the 
conservative approach would be to use a smaller fmax to 
stay away from non-linear scales. However, we feel this is 
too conservative as the ratio R is f-independent and the 
power spectra that enter the error calculation are known 
with quite good accuracy that will improve considerably 
with future observations. 



FIG. 4: Dark energy Icr constraints on wo and Ude using 
the same assumptions as in Fig.j^ We show a comparison of 
the LSST-only results (dashed blue contour, labelled “LSST”) 
with the results of combining SKAl-Mid for the foreground 
lenses and LSST for the background sources (dotted black 
contour, labelled “SKAl(j')”) using ^max = 10"^. Results using 
SKA2-Mid for the foreground lenses are also shown (solid red 
contour, labelled “SKA2(^)”). 

The Icr marginalised results for the constant equation 
of state parameter uncertainty are 

- <5u;oltT=30oo - 0.12, 

^-olL^:;LTo^-^tuo|LT=io^-0.10, 

(14) 

Here we note that our results for the LSST-only case with 
•^max = 10 "*^ are in agreement with previously published 
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geometrical constraints, although an exact comparison is 
difhcult to make as different works have used different 
methods, assumptions, approximations and priors. Our 
approach is closer to the one followed by [3]. There, a 
constraint Swq ^ 0.1 is found for a survey with parame¬ 
ters similar to LSST. 

Using SKAl for the foreground measurement offers 
comparable constraints to using LSST, but the cross¬ 
correlations could be less susceptible to systematic con¬ 
tamination; improving the observations assuming SKA2 
parameters can improve the constraints. However, the 
improvement is not very large as the dominant terms in 
the Fisher matrix calculations are the LSST shape noise 
terms shown in Fig. bottom panel. 


C. Results using a LSST-like survey for (/, 6i) and 
a SKA-Mid-like Intensity Mapping survey for (62) 

We next consider the case where we probe the fore¬ 
ground (/) lenses and the (&i) sources with a galaxy 
survey like LSST, but we use a SKA-Mid-like Intensity 
Mapping survey for the (& 2 ) sources. In order to calcu¬ 
late the noise in the k measurement for this case, we use 
the quadratic estimator technique developed in [ISIIIII- 
This estimator takes into account the discreteness and 
clustering of galaxies and is applicable to any redshift be¬ 
low reionization. Using the intensity mapping approach 
and an SKA-like array, the method performs better than 
21cm galaxy surveys at high redshifts. Hence, in the fol¬ 
lowing we are going to apply the weak leasing intensity 
mapping estimator for the second background sources 
( 62 ) at various redshifts z > 2, but still use LSST for 
the closer background sources ( 61 ) at z ~ 1. We are also 
going to use LSST for the foreground lenses (/), but keep¬ 
ing in mind that an SKA-like interferometer can perform 
equally well or better, as shown in Fig. 

For a detailed discussion and derivation of the results 
used here, see |27j . The expression for the leasing re¬ 
construction noise Nb^ {£) from the intensity mapping 
method is rather lengthy, so we will not include it here, 
but the interested reader is referred to [ 22 ], Appendix C. 
To summarise, Nb^ involves the underlying dark matter 
power spectrum Pss, the HI density Uhi(z:) as well as the 
HI mass (or luminosity) moments up to 4th order and, of 
course, the thermal noise of the instrument (in our case 
the SKA). In ^7} it was found that the signal-to-noise 
is strongly dependent on the possible evolution of the 
HI mass function. More specifically, assuming the no¬ 
evolution scenario (which is the most conservative, but 
also less realistic approach), precise measurements can 
be made with an SKA2-like instrument; however if we 
assume instead a model where the HI density Uhi(z) in¬ 
creases by a factor of 5 by redshift z = 3 and then slowly 
decreases towards redshift z = 5, as suggested by the 
DLA observations from [3S], high signal-to-noise can be 
achieved even with SKAl. 

The parameters for SKA2-Mid are the same as in the 


previous subsection. The source redshift is Zbj = 3 with a 
bandwidth of 40 MHz. The results are shown in Figure 
and we have used ^max = 10^- The largest contour corre¬ 
sponds to the no-evolution HI scenario, while the smaller 
one assumes the aforementioned HI evolution model. 



FIG. 5: Dark energy Icr constraints on wo and Ude Dark en¬ 
ergy la constraints on wq and Ddb using the same assump¬ 
tions as in Fig.|^ We have used a combination of a LSST-like 
galaxy survey in the foreground and the 61 background, and 
a SKA2-Mid-like IM survey in the 62 background with pa¬ 
rameters described in the main text and ^max = 10^. Sources 
are at zj^ = 3. We show results assuming no HI evolution 
(larger contour) and using the evolution scenario described in 
the main text (smaller contour). 

The results are again competitive with the LSST-only 
case, and they considerably benefit from a possible HI 
evolution. The result for HI sources at z = 3 using the 
evolution model and SKA2-Mid is Swq ~ 0.08. 

D. Results using a LSST-like survey for (/, for) and 
21 cm lensing from the Epoch of Reionization (62) 

Let us now consider the case where our second back¬ 
ground ( 62 ) sources are at the Epoch of Reionization 
(EoR). This provides a longer lever arm for the ratio 
R and uses information from an era that will be ex¬ 
plored for the first time with the SKA. The possibility 
of measuring the lensing signal with 21 cm emission from 
the EoR has been studied previously. In and [30] 
the convergence estimator and the corresponding lensing 
reconstruction noise were calculated assuming that the 
temperature (brightness) distribution is Gaussian, which 
is a reasonable approximation at the EoR, at least while 
the ionised regions are small. An important advantage of 
21 cm lensing is that one is able to combine information 
from multiple redshift slices. 

In Eourier space, the temperature fluctuations are di¬ 
vided into wave vectors perpendicular to the line of sight 
kj^ = l/D, with V the angular diameter distance to the 
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source redshift, and a discretised version of the paral¬ 
lel wave vector /cy = where £ is the depth of the 
observed volume. Considering modes with different j 
as independent, an optimal estimator can be found by 
combining the individual estimators for different j modes 
without mixing them. The three-dimensional leasing re¬ 
construction noise is then found to be [53] 


= (£ 4 / 4 ) X 

r dH' [I'-lQ^.+l- 

UJ (2^)" 


where 

and 

I /^N 

"T • 


(15) 

(16) 

(17) 


Following this approach, we consider the case where 
we probe the foreground (/) lenses and the ( 6 i) sources 
with a galaxy survey like LSST, and we use a SKA-like 
instrument for the ( 62 ) sources at the EoR redshift, which 
we assume to be instantaneous at order to 

facilitate the comparison with the previous Sections, we 
are initially going to consider a SKA2-Mid-like instru¬ 
ment but with a smaller observation bandwidth equal to 
8 MHz. That is because the Gaussian quadratic leasing 
estimator is optimised for the case where the statistical 
properties of the 21 cm radiation signal and noise are con¬ 
stant within a band, and an observation bandwidth of a 
few MHz is small enough so that this assumption is jus¬ 
tified. The results are shown in Figure for £niax = lO^ 
(solid black contour). They give a marginalised uncer¬ 
tainty <5u>o — 0 . 11 . 

However, SKA-Mid is not designed to probe low 
enough frequencies (equivalently, high redshifts), and the 
current plans for SKA-Low which will observe the EoR 
redshifts will scan a very small sky area hence it will not 
be able to achieve the required precision in the k measure¬ 
ment. In dillT] it was found that a compact SKA-like 
instrument, i.e. an array with smaller maximum baseline 
and a smaller core collecting area, would perform a weak 
lensing intensity mapping survey equally well or even bet¬ 
ter than the current SKA-Mid design. This means that 
the precision required could be achieved by using a more 
compact configuration. We therefore present results con¬ 
sidering a SKA-like “purpose-built” compact interferom¬ 
eter. The collecting area is Acou = 0.5 km^ and the max¬ 
imum baseline Utei = 4 km, keeping the rest of the pa¬ 
rameters the same. The results are shown in Figure 
(dashed red contour). With a more compact configu¬ 
ration, we can obtain identical constraints from cross¬ 
correlation cosmography with significantly less collecting 
area (Aeon = 0.5 km^, compared to Acoii = 0.8 km^ as¬ 
sumed in the original SKA2 configuration.) 



FIG. 6 : Dark energy Icr constraints on wo and Hde using 
the same assumptions as in Fig. We have used a combi¬ 
nation of a LSST-like galaxy survey in the foreground and 
the 61 background, and a SKA2-Mid-like instrument in the £2 
background at the EoR redshift assumed to be z = 7 (solid 
black contour). We also show results with a more compact 
SKA-like instrument (dashed red contour). 


E. Results using a LSST-like survey for (/, foi) and 
CMB lensing with a COrE-like satellite for (62) 


Finally, we consider the case where we use the CMB as 
the second source plane. This possibility has been stud¬ 
ied in the past EH 131] and its main advantages are the 
very well determined source redshift and distance and the 
fact that it provides the largest possible lever arm for the 
distance ratio. For the CMB case, the two-dimensional 
lensing reconstruction noise is found to be [33] 


iV(£)|cmb = (£V4)X 


d^e [i'-icTT + i-(i-i')qy_J 


(27r)2 


O/^tot ^tot 


(18) 


where is the CMB TT power spectrum and 


cr = cr+ ( 19 ) 

with Cf the instrumental noise for the CMB survey. We 
will consider a future COrE-like satellite with FWHM 
a = 3.0' and temperature noise Ay = 1 /rK'. Then we 
have 


exp[£(£ -b l)(TV81n2]. ( 20 ) 

In Figure [^ we show the derived constraints using the 
aforementioned satellite for the second source plane ( 62 ) 
and the LSST-like survey for (/,5i). The marginalised 
wq uncertainty is found to be 5wo ~ 0.22. Comparing 
these results to the SKA-like EoR case from the previ¬ 
ous Section we see that the latter give much better con¬ 
straints. That is because the CMB reconstruction comes 














with higher statistical errors compared to the 21cm re¬ 
construction, and the leverage we get from the higher 
redshift is not enough to compensate. 



FIG. 7: Dark energy la constraints on wo and Ude using 
the same assumptions as in Fig. We have used a combi¬ 
nation of a LSST-like galaxy survey in the foreground and 
the for background, and a COrE-like CMB survey in the 62 
background. 


F. Noise terms comparison 

To consolidate our results, we compare the noise terms 
used in the Fisher matrix calculations for the various sur¬ 
vey combinations we have considered. The top panel of 
Fig. [8] compares the tracer density power spectra and 
noise terms for LSST (dashed green lines), SKAl (solid 
black lines) and SKA2 (dotted red line) at the 2 ;/ = 0.5 
foreground slice. Cs^Sg and Cs^^Shi are comparable but 
not exactly the same because their redshift binnings 
differ slightly; for LSST, Az = 0.1 but the choice of 
A/ = 40MHz bandwidth for the intensity mapping sur¬ 
vey corresponds to a narrower bin at Zf = 0.5. The 
LSST and SKAl noise terms (JV/) are comparable, while 
the SKA2 noise is much lower. This explains the results 
of Fig. where we found comparable constraints from 
LSST-only and SKAl(f) cases, and stronger constraints 
using SKA2 to map the foreground density. 

For the intermediate lensing data, which we typically 
assume is centred at z/,^ = 1.3 to maximise the distance 
to the higher redshift lensing, the HI lensing noise is 
not competitive with that of the LSST. However, the 
HI intensity mapping measurements can be competitive 
for highest redshift sample, where they can observe at 
higher redshifts than are possible with LSST. The middle 
panel of Fig. compares the convergence power spectra 
and Nb^ noise terms for LSST at Zb^ = 2.4 (solid black 
lines), SKA2 at Zb^ = 3 assuming no HI evolution (dashed 
blue lines) and SKA2 at Zb 2 = 3 assuming the HI evo¬ 


lution model described in the main text (dotted-dashed 
red line). As we can see, the LSST and SKA2 HI evolu¬ 
tion noise terms are comparable, and in the SKA2 case 
we also have a longer lever arm because the &2 sources 
are at higher redshift. 





FIG. 8 : A comparison of the noise terms for the various survey 
combinations considered. See text for further details. 

Finally, the bottom panel of Fig. [^compares the con¬ 
vergence power spectra and noise terms for a SKA2- 
like instrument at the EoR redshift assumed to be 2 : = 7 
(solid black lines) and a COrE-like satellite with sources 
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at the CMB redshift z = 1090 (dashed magenta lines). 
We see that the noise level of the COrE-like mission is 
higher than SKA2 and quickly diverges as we reach the 
limits (£max — 3600), set by the beam resolution. Again, 
this is a consequence of the CMB providing only a sin¬ 
gle source plane, while the HI data offers multiple source 
planes that can be co-added to reduce the noise. 

IV. CONCLUSIONS 

In this paper, we have shown how the HI intensity map¬ 
ping technique can be used to enhance cross-correlation 
cosmography studies. Cross-correlation cosmography is 
based on the scaling of the cross-correlation signal 
with the redshift of the background source population. 
The ratio of the signals from the same foreground lens 
population to two different background populations sim¬ 
plifies to a geometrical distance ratio that only depends 
on dark energy parameters and curvature. Intensity map¬ 
ping is a technique that treats the 2 Icm emission as a 
continuous unresolved background, without resolving or 
even identifying (in angular resolution, not frequency) 
individual galaxies. It offers excellent redshift resolution 
and a longer lever arm for the background sources. 

We performed a comprehensive study of the possibility 
of combining optical galaxy surveys, in particular LSST, 
with HI intensity mapping surveys using a SKA-like in¬ 
strument in order to derive geometrical dark energy con¬ 
straints. Our results show that using the SKA to measure 
the foreground density field / and/or the background 
convergence field 62 at high redshifts z > 2 can be benefi¬ 
cial. More specifically, we find that a constant equation of 
state for dark energy can be constrained to ~ 8 % for a sky 
coverage /sky = 0.5 and assuming a ct(Hde) = 0.03 prior 
for the dark energy density parameter. One major uncer¬ 
tainty is the unknown evolution of the HI density param¬ 
eter and the form of the HI mass function that is crucial 
for the modelling of the lensing reconstruction noise us¬ 
ing intensity mapping. However, the no-evolution model 
we considered is the most conservative scenario. Using 
a more optimistic - but also more realistic - evolution 
scenario significantly improves the expected constraints. 
The Epoch of Reionization itself can potentially be used 
as our second background source; it provides a longer 
lever arm for the distance ratio and the combination of 
the contribution from many redshifts slices results in a 
low noise level, offering significantly better constraints 
compared to using the CMB, even when assuming COrE- 
like satellite observations. 

Looking at Table [l] which summarizes all the survey 
combinations we have considered, one notices that we 
never consider all (/, 61 , 62 ) populations to be probed 
with intensity mapping; we always use an LSST-like sur¬ 
vey for the closer background population hi. That is 
because the weak lensing intensity mapping estimator is 


competitive with optical galaxy surveys at high redshifts 
z > 2 , where we also have the advantage of a longer lever 
arm. In lower redshifts a galaxy survey performs better. 
We also note that in cases HI.C, HI.D and HI.E in Ta- 
blej^we can replace LSST with SKAI for the foreground 
lenses (/), as they perform equally well (see Fig. |^. 

Foreground contamination can be a significant con¬ 
cern. Foregrounds are the most important source of 
systematic error when it comes to the intensity map¬ 
ping technique, but interferometric techniques can filter 
them out. In particular for our lensing studies, where the 
frequency-dependent foreground contribution is large, it 
has been shown that the foreground subtraction tech¬ 
niques will remove parallel k modes |23L I29j . meaning 
that the lensing reconstruction noise will increase some¬ 
what. However, if the subtraction techniques are suc¬ 
cessful, only the first few modes will be removed and the 
signal-to-noise of the measurements will not be signifi¬ 
cantly affected m- 

To conclude, we have demonstrated that HI inten¬ 
sity mapping can be used for cosmographic studies, to 
complement and compete with the state-of-the art op¬ 
tical galaxy surveys. In addition, the weak lensing in¬ 
tensity mapping technique performs very well across a 
wide range of post-reionization redshifts, and by using 
tomographic information in the measured convergence 
it is possible to infer how the matter power spectrum 
and the growth function evolve with time [57]. These 
two avenues (cosmography and lensing tomography) can 
be explored simultaneously in order to derive combined 
growth-geometry constraints on the dark energy param¬ 
eters. An example of how this can be done is presented 
in [3|. The dark energy parameters can be split into two 
kinds, those that enter the growth factor, and those that 
enter the geometrical distances. Such a splitting allows a 
robust consistency check if the equation of state w values 
obtained separately from geometry and from growth are 
in agreement; if, however, they disagree, we can identify 
and remove important systematics or contaminations, or 
revise incorrect assumptions about the behaviour of the 
mass fluctuations. The results of such a consistency test 
are shown in |3]. 

Our studies would benefit from the modelling of a 
“realistic”, non-instantaneous reionization history, from 
more accurate measurements of Hhi6hi as a function of 
redshift and from detailed simulations of foreground con¬ 
tamination and subtraction techniques. These investiga¬ 
tions are under way and will be presented in future work. 
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